Infrared phonons and specific heat in the gapped quantum magnet Ba 3 Cr 2 08 
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We report on the phonon spectrum of the spin-gapped quantum magnet BaaC^Os determined 
by infrared spectroscopy, and on specific heat measurements across the Jahn- Teller transition in 
magnetic fields up to 9 T. Phonon modes split below the Jahn- Teller transition, which occurs at 
Tjt — 70 K as detected by specific-heat measurements. The field-dependent specific heat data is 
analyzed in terms of the contributions from lattice, magnetic and orbital degrees of freedom. In 
contrast to the isostructural compound S^C^Os our analysis does not indicate the existence of 
orbital fluctuations below the Jahn- Teller transition in BasC^Os. 

PACS numbers: 78.30.-j,71.70.-d,75.10.Hk 
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I. INTRODUCTION 

Spin-gapped quantum antifcrromagncts with dimer- 
ized ground states have been intensively studied due to 
the observation of various exotic phenomena in these 
systems^— Based on dimers of Cr 5+ (3d 1 , 5 = 1/2) or 
Mn 5+ (3d 2 , 5 = 1) ions, a series of isostructural dimer- 
ized compounds, such as Ba3Cr 2 08r r — Sr 3 Cr 2 08^ and 
Ba 3 Mn 2 0g^ exhibits a magnon-condensation state with 
the application of an external magnetic fields 

The room-temperature crystalline structure of these 
compounds is hexagonal with space group R2>m^ i 10 : 11 In 
Ba3Cr208 and S^C^Os, a tetrahedral crystal field splits 
the 3d levels of the Cr ions into lower-lying doubly de- 
generated e orbitals and higher-lying triply degenerated 
t orbitals. Thus, the systems are Jahn- Teller active and 
they are expected to undergo a structural phase transi- 
tion to remove the orbital degeneracy. Neutron powder 
diffraction indeed revealed such a structural transition 
at 70 K in Ba 3 Cr 2 8 and at 275 K in Sr 3 Cr 2 8 from 
hexagonal to monoclinic with space group C2/ 

The magnetic susceptibility of Ba3Cr 2 0s has been re- 
ported by Nakajima et al.r^ and Aczcl et alJ^- They 
have shown that the susceptibility can be described 
by a Bleaney-Bowers modelj 1 ^ where they considered a 
magnetic structure consisting of weakly-interacting spin 
dimers of Cr 5+ ions. Fitting to this model resulted 
in an intra-dimer exchange interaction Jo = 2.15 meV 
consistent with the value of 2.38 meV determined by 
neutron experiments^ The detailed magnetic structure 
has been studied theoretically based on extended Hiickel 
tight binding calculations The authors of Ref. UH have 
examined all the possible exchange paths between the 
nearest-neighboring Cr ions. They found that the ex- 
change interaction between adjacent Cr 5+ ions in the 
hexagonal c/j direction is most important, while the other 
exchange interactions are much smaller. This was con- 
firmed by inelastic neutron experiments in Ba3Cr 2 0s^ 



The magnetic structure of Sr3Cr 2 0s is the same as 
Ba3Cr 2 0g, but the intradimer exchange interaction of 
Sr3Cr 2 C>8 is 5.55 meVj 1 ^ about twice the corresponding 
value in Ba3Cr 2 0s, and the Jahn- Teller transition tem- 
perature 285 K in Sr3Cr 2 0§ is about four times higher 
than that in Ba3Cr 2 0gri^ In addition, an anomalous 
strong damping of the optical phonon modes has been ob- 
served down to about 120 K in Sr3Cr 2 0s, which is below 
the Jahn- Teller transition temperature of Sr3Cr 2 0g i 17 i 18 
Moreover, an analysis of the specific heat revealed a resid- 
ual orbital entropy below the structural transition, indi- 
cating strong fluctuations in the orbital degrees of free- 
dom. 

In this work we investigated the polar phonons and the 
specific heat in Ba3Cr 2 0s to search for a similar fluctua- 
tion regime. We have resolved the phonon spectra from 
the reflectivity measurements and found the splitting of 
phonon modes right below the Jahn- Teller transition, in- 
dicating that strong fluctuations of the lattice degrees 
of freedom are absent in Ba3Cr 2 0s- This clearly distin- 
guishes the system from Sr 3 Cr 2 0s. The analysis of the 
specific heat confirms this difference between the com- 
pounds, because the orbital entropy associated with the 
splitting of the orbital ground state is completely recov- 
ered at the Jahn- Teller transition in Ba3Cr 2 0s. 



II. EXPERIMENTAL DETAILS 

Single crystals of Ba 3 Cr 2 0s were grown by the 
floating-zone method. They were cleaved perpendic- 
ular to the hexagonal c/j-axis and checked by single- 
crystal x-ray diffraction. The heat capacity was mea- 
sured in a Quantum Design physical properties measure- 
ment system at temperatures 1.8 < T < 300 K with 
applied magnetic fields up to 9 T. Reflectivity measure- 
ments were carried out on optically polished samples for 
4 < T < 300 K in the far- and mid-infrared range using 
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the Brukcr Fourier-transform infrared spectrometers IFS 
113v and IFS 66v/S with a He-flow cryostat (Cryovac). 



III. EXPERIMENTAL RESULTS AND 
DISCUSSION 



lowest lying mode E U {1) could not be completely resolved 
and only the high frequency flank has been fitted. Thus, 
the fitting parameters for E u (l) contain a larger experi- 
mental uncertainty. The fitting parameters for the 7E U 
modes are listed in Table HI 



A. Infrared spectra 

Figure [T] shows reflectivity spectra of BaaC^Os mea- 
sured for the radiation electric field E parallel with the 
crystalline hexagonal a^-planc (E |j a/j&^-plane) at 
295 K and 4 K, which are above and below the Jahn- 
Teller transition temperature Tjt = 70 K, respectively. 
The seven modes observed above Tjt are tentatively as- 
signed as E u (j) (j=l,.,.,7). The primitive unit cell of the 
high-temperature hexagonal structure with space group 
R3m (No. 166) 19 i 20 contains one chemical formula unit 
and, hence, a total of 39 normal modes exist. The IR 
active modes are characterized by the irreducible rep- 
resentations 6A 2 „(E || c h ) + 7E U [E || (a h ,b h )]&2L The 
number of normal modes is in agreement with the obser- 
vation of infrared and Raman active phonons reported 
for the isostructural system SrgC^Osi 17 ' 18 The overall 
appearance of the BaaCi^Os far-infrared spectra are sim- 
ilar to the ones of SraC^Os, but the BaaC^Os phonons 
are shifted to lower cigcnfrcqucncies consistent with the 
larger molar mass of BaaC^Os^ The absolute value of 
the reflectivity is relatively low, which might be due to 
the imperfection of the sample surface, and similar val- 
ues have been found for S^C^Os-— The weak features 
in the frequency range between modes E u (6) and E u {7) 
are also assigned to originate from defects on the sample's 
surface. 

The quantitative analysis of the phonon modes is 
performed by fitting the reflectivity spectrum with the 
Drude-Lorentz model. According to this model, the 
phonon contribution of the complex dielectric function 
is given by 
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where each phonon mode j is described by the eigen- 
frequency u>oj, damping coefficient jj and effective ionic 
plasma frequency flj. The parameter takes into ac- 
count contributions of higher-lying electronic excitations 
and was used as a free fitting parameter ranging between 
~ 4 at room temperature and ~ 7 just above the Jahn- 
Teller transition. At normal incidence, the reflectivity 
R(u>) is related to e(u>) via 



R(cj) 



(2) 



The fit result for the room-temperature spectrum using 
seven oscillators is given by the solid line in Fig. [TJ which 
agrees well with the experimental data. Note that the 



TABLE I: Eigenfrequency uq, effective ionic plasma frequency 
Q. and damping coefficient 7 of all 7E U phonon modes in 
BaaCraOg at 295 K. 



E u 




7 Phonons 




Mode 


wo (cm x ) 


n (cm -1 ) 


7(cm i ) 


1 


52.6 


101.2 


13.5 


2 


119.7 


187.9 


14.0 


3 


174.7 


182.8 


10.6 


4 


188.4 


174.1 


26.2 


5 


274.7 


166.0 


172.1 


6 


330.5 


120.4 


14.9 


7 


759.2 


718.2 


45.5 



For the monoclinic C2/c (No. 15b) structure below the 
JT transition, the number of expected normal IR modes 
increases to 19A U for the radiation electric field parallel 
with monoclinic b m axis and 20B U for the electric field 
parallel with monoclinic a m c m -plane, i.e. T = 19^4 U (E | 
b m ) + 20£> U [E || (a m , c m )].— >2i According to the geomet- 
ric relation between hexagonal and monoclinic axes ah = 



a m -b m ), b h = -h(a m + b m ), and c h = \c m -\a m de- 



termined by neutron diffraction measurements 1 1 1 1 16 the 
IR spectra measured with E || a^6h-plane below Tjt 
should display all the 19^4 U and 20S U modes. However, 
the coexistence of three monoclinic twins in the single 
crystal^i^ and the low monoclinic symmetry make it dif- 
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FIG. 1: (Color online) Reflectivity spectra of BasC^Os mea- 
sured with E || dhfrh-plane at 295 K and 4 K. The spectrum 
of 4 K is shifted with respect to that of 295 K by a constant 
of 0.1. Infrared-active phonon mode E u (j)(j = 1, 7) are 
marked at the spectrum of 295 K for the hexagonal structure. 
Solid line is a fit of Drude-Lorentz model to the data of 295 K 
as described in the text. 
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FIG. 2: (Color online) Reflectivity spectra corresponding to 
IR phonon modes E u (2) and E u {7) measured with E || anbh- 
plane at various temperatures. The spectra are shifted with 
respect to that at 295 K by a constant in order to clearly il- 
lustrate the evolution of phonon modes with decreasing tem- 
perature. 



ficult to determine the eigcnfrcqucncy, damping coeffi- 
cient and plasma frequency of the phonon modes below 
the Jahn- Teller transition^ 

The temperature-dependent reflectivity spectra reveal 
that the hexagonal mode of E u (2) and E u (7) split below 
Tjt, while a splitting of the other modes, which arc hard- 
ening with decreasing temperature, could not be resolved 
below Tjt- Therefore, we focus on the temperature- 
dependent features of modes E u {2) and E u {7), and try 
to describe the low-temperature modes with Lorentzian 
functions according to Eqs. (JTJ and (|2J). As indicated by 
the arrows in Fig. we model the low-temperature spec- 
tra with two oscillators for both the E u (2) and the E u {7) 
modes. For the monoclinic structure, the phonon modes 
cannot be uniquely determined in a single polarization 
measurement. Therefore, the obtained fit parameters in- 
cluding do not necessarily reflect the intrinsic dielec- 
tric properties of the system and it might be misleading 
to show the dielectric parameters. Instead, the reflectiv- 
ity spectra should be measured with three different polar- 
izations so as to determine the three independent matrix 
elements of the dielectric matrix for E || a m c m -planei22 
This is infeasible for BaaCr^Os with the coexistencce of 
three monoclinic twins^^ Thus the modeling here with 
Lorentzian functions is rather a way to parametrize the 
reflectivity spectra and reveal clear changes occurring at 
Tjt- 

The fitting results of eigenfrequency, damping coeffi- 
cient and plasma frequency are summarized as a func- 
tion of temperature in Fig. [3l which clearly shows the 
splitting of phonons right below Tjt, consistent with the 
phase transition to monoclinic structure. This is an ev- 
ident contrast to the case of S^C^Os, which is an iso- 
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FIG. 3: (Color online) Temperature dependence of eigenfre- 
quency (a, d), damping coefficient (b, e) and ionic plasma 
frequency (c, f) for modes E u (2) and E u (7) on a semi- 
logarithmic scale. The Jahn- Teller transition temperature 
Tjt is indicated by vertical dashed line. Open symbols signify 
the results of parametrization below Tjt (see text). 

structural compound of BaaC^Os^ In Sr 3 Cr 2 08, the 
new phonons are not emerging right below the Jahn- 
Teller transition due to strong fluctuations. With de- 
creasing temperature, the fluctuations are reduced and 
all the modes of monoclinic structure are observed when 
the temperature is lower than 100 KJ^ 

The modes marked by "+" in Fig. [5] at 812 cm -1 and 
834 cm -1 are visible right at 40 K and lower tempera- 
tures. These modes could be new phonon modes of the 
monoclinic structure, which have small spectral weight 
and large damping effect that could not be distinguished 
from the strong modes corresponding to E u (7) at higher 
temperatures. 



B. Specific heat 

Figure HJa) shows the specific heat of BaaCr^Os mea- 
sured as a function of temperature. Two anomalies can 
be observed in the specific-heat data. A sharp anomaly 
indicated by a dashed vertical line can be easily recog- 
nized at 70 K. Another feature observed below 20 K is 
marked by an arrow. The inset of Fig. HJa) shows the 
specific heat over temperature C/T vs. T measured in 
various external magnetic fields. The anomaly at Tjt = 
70 K is independent on the magnitude of magnetic field, 
and marks the structural phase transition from hexagonal 
R3m to monoclinic C2/c, in agreement with the results 
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FIG. 4: (Color online) (a) Specific heat of BaaC^Os without 
applied magnetic field. Solid line is a fit of lattice contribution 
to specific heat according to Debye and Einstein models as 
described in the text. Inset shows specific heat divided by 
temperature C/T measured at various magnetic fields, (b) 
Residual specific heat (C — Ciatt) /T and (c) the corresponding 
entropy obtained at different fields, (d) Magnetic contribution 
to specific heat fitted by considering the Zeeman interaction 
of isolated spin dimers with magnetic field (solid lines) as 
described in the text. 



obtained from neutron diffraction measurements.— The 
anomaly below 20 K is very sensitive to the magnitude 
of the magnetic field indicating that is related to the spin 
degrees of freedom. 

The lattice contribution to the specific heat Ci a tt can 
be described according to the phononic spectra that have 
been determined in Section IIII Al We model the lattice 
contribution by a sum of one isotropic Debye (D) and 
four isotropic Einstein terms (£1.2,3.4) with the fixed ra- 
tio D : E 1 : E 2 : E 3 : £4 = 1 ': 3 : 4 : 3 : 2 be- 
tween these terms accounting for the 39 degrees of free- 
dom per formula unit. This ratio has been previously 
used to describe the specific- heat data of Sr 3 Cr 2 08.— 
The resulting contribution to the specific heat shown as 
a solid line in Fig. @|a) has been obtained with the Debye 



and Einstein temperatures Od = 123 K, B E \ = 137 K, 
6 E2 = 265 K, 9 E3 = 658 K, and 6 Ei = 1270 K in 
agreement with the frequency ranges where 1R active 
phonons of the hexagonal structure occur. The fitting 
curve agrees well with the specific-heat data above Tjt, 
while below Tjt it deviates from the experimental data 
indicating that the magnetic and orbital contributions 
to the specific heat are important below the Jahn- Teller 
transition. Since BasC^Os has larger molecular mass 
than Sr 3 Cr 2 8 , the obtained Debye and Einstein tem- 
peratures of Ba 3 Cr 2 8 should be lower than those of 
Sr3Cr20g The fitted Debye-temperature ratio between 
Ba3Cr20g and S^C^Os is 0.91, which is consistent 
with the ratio yfmsco ■ y/rn bco = 0.88, where msco 
and tubco are the molecular masses of S^C^Os and 
BaaC^Os, respectively. 

The residual specific heat divided by temperature (C — 
Ciatt) /T is shown in Fig.@|b). One can clearly recognize 
a A-shaped anomaly at the Jahn- Teller transition around 
70 K and the magnetic anomaly below 20 K, which are 
well separated in the temperature scale. The correspond- 
ing entropy S(T) = f Q T dd{C — Ciatt)/® reaches a value 
slightly lower than the expected S = 4i?ln2 for the 
sum of magnetic and orbital contribution as shown in 
Fig.[4|c), where R is the molar gas constant. For a spin- 
1/2 dimer, the entropy increases by 2 In 2 from the or- 
dered to the disordered spin state. Because there is essen- 
tially one spin dimer per formula unit of BasC^Os, the 
magnetic contribution to the entropy should be 2i?ln2. 
This is reached at 20 K indicated in Fig. @Jc) . At this 
temperature, the specific heat measured at different fields 
begins to deviate from each other on cooling. Above the 
Jahn- Teller transition the e orbitals of each Cr 5+ ion are 
degenerate and the orbital degrees of freedom will con- 
tribute to the entropy by another 2i?ln2, where the two 
Cr ions per formula unit have been taken into account. 

In order to describe the magnetic contribution to the 
specific heat, we consider a model with the intra-dimer 
exchange interaction and the Zeeman interaction but 
ignore the inter-dimer exchange interactions which are 
much smaller than the intra-dimer interaction^ Thus, 
we can model the magnetic contribution to the specific 



heat by C mag (T) = using E 
the partition function Z = 



1 V 3 

Z 



-P ei with 



-fa 



and the energies 



£0,1,2,3 = 0, Jo - 9(J.bH, J , J + gfJsH, where [i B is the 
Bohr magneton, H is the magnetic field, and (3 = l/ksT 
with the Boltzmann constant fc^. The solid lines in 
Fig. Hid) are the calculated specific heat curves according 
to this model, where the intra-dimer exchange interac- 
tion Jo = 2.38 meV was determined by inelastic neutron 
experiments^ and the g-factor g = 1.94 was measured by 
electron spin resonance experiments^ This simple model 
captures the most important features of the experimental 
data. The modeled specific heat reaches the maximum 
nearly at the same temperature as that of the experi- 
mental data. The maximum shifts to lower temperature 
and becomes smaller with the increase of magnetic field, 
since the gap between the lower-lying triplet state and the 
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singlet state is reduced due to Zeeman interaction with 
the magnetic field. Because inter-dimer interactions are 
not considered in the model, the modeled specific heat 
should be smaller than the experimental data as observed 
in Fig. HJd). Note that this discrepancy could also orig- 
inate to a certain amount from the uncertainties of the 
changes in the phonon density of states in the low tem- 
perature structure. 

IV. CONCLUSION 

In summary, the infrared phonon spectrum measured 
with the polarization perpendicular to the hexagonal 
c^-axis in BaaC^Os exhibits all expected 7E U modes 
above the structural phase transition at 70 K. Below this 
transition, which is associated with a Jahn- Teller dis- 
tortion and orbital ordering, a splitting of some phonon 
modes and the appearance of new modes indicate the re- 
duced symmetry in the distorted phase. In the specific 
heat a clear anomaly is visible at 70 K, which does not 
change with applied magnetic field. An additional low- 
temperature anomaly in the specific heat is assigned to 
originate from magnetic excitations and its dependence 
on magnetic field can be explained by considering the 
thermal population of triplet states and their Zeeman 



splitting in external magnetic fields. The entropy change 
associated with orbital ordering is found to be reached 
just above the structural Jahn- Teller transition without 
any sign of orbital fluctuations below 70 K. In compar- 
ison with the isostructural compound S^C^Os, where 
the polar phonons and the specific heat bear signatures of 
strong fluctuations and competing interactions, the fact 
that this is not the case for BaaC^Os might help to 
understand which interactions are responsible for the cu- 
rious behavior of Sr 3 Cr 2 08 and single out the differences 
in these otherwise very similar systems. 
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